Systems Biology Graphical Notation: Entity Relationship language Level 1

Standard graphical representations have played a crucial role in science and engineering throughout the last century. Without electrical symbolism, it is very likely that our industrial society would not have evolved at the same pace. Similarly, specialised notations such as the Feynmann notation or the process flow diagrams did a lot for the adoption of concepts in their own fields. With the advent of Systems Biology, and more recently of Synthetic Biology, the need for precise and unambiguous descriptions of biochemical interactions has become more pressing. While some ideas have been advanced over the last decade, with a few detailed proposals, no actual community standard has emerged. The Systems Biology Graphical Notation (SBGN) is a graphical representation crafted over several years by a community of biochemists, modellers and computer scientists. Three orthogonal and complementary languages have been created, the Process Descriptions, the Entity Relationships and the Activity Flows. Using these three idioms a scientist can represent any network of biochemical interactions, which can then be interpreted in an unambiguous way. The set of symbols used is limited, and the grammar quite simple, to allow its usage in textbooks and its teaching directly in high schools. The first level of the SBGN Entity Relationship language has been publicly released. Shared by the communities of biochemists, genomicians, theoreticians and computational biologists, SBGN languages will foster efficient storage, exchange and reuse of information on signalling pathways, metabolic networks and gene regulatory maps

Model storage, exchange and integration

The field of Computational Systems Neurobiology is maturing quickly. If one wants it to fulfil its central role in the new Integrative Neurobiology, the reuse of quantitative models needs to be facilitated. The community has to develop standards and guidelines in order to maximise the diffusion of its scientific production, but also to render it more trustworthy. In the recent years, various projects tackled the problems of the syntax and semantics of quantitative models. More recently the international initiative BioModels.net launched three projects: (1) MIRIAM is a standard to curate and annotate models, in order to facilitate their reuse. (2) The Systems Biology Ontology is a set of controlled vocabularies aimed to be used in conjunction with models, in order to characterise their components. (3) BioModels Database is a resource that allows biologists to store, search and retrieve published mathematical models of biological interests. We expect that those resources, together with the use of formal languages such as SBML, will support the fruitful exchange and reuse of quantitative models

JSBML: a flexible Java library for working with SBML

The specifications of the Systems Biology Markup Language (SBML) define standards for storing and exchanging computer models of biological processes in text files. In order to perform model simulations, graphical visualizations and other software manipulations, an in-memory representation of SBML is required. We developed JSBML for this purpose. In contrast to prior implementations of SBML APIs, JSBML has been designed from the ground up for the Java™ programming language, and can therefore be used on all platforms supported by a Java Runtime Environment. This offers important benefits for Java users, including the ability to distribute software as Java Web Start applications. JSBML supports all SBML Levels and Versions through Level 3 Version 1, and we have strived to maintain the highest possible degree of compatibility with the popular library libSBML. JSBML also supports modules that can facilitate the development of plugins for end user applications, as well as ease migration from a libSBML-based backend

Structural Analysis and Stochastic Modelling Suggest a Mechanism for Calmodulin Trapping by CaMKII

Activation of CaMKII by calmodulin and the subsequent maintenance of constitutive activity through autophosphorylation at threonine residue 286 (Thr286) are thought to play a major role in synaptic plasticity. One of the effects of autophosphorylation at Thr286 is to increase the apparent affinity of CaMKII for calmodulin, a phenomenon known as “calmodulin trapping”. It has previously been suggested that two binding sites for calmodulin exist on CaMKII, with high and low affinities, respectively. We built structural models of calmodulin bound to both of these sites. Molecular dynamics simulation showed that while binding of calmodulin to the supposed low-affinity binding site on CaMKII is compatible with closing (and hence, inactivation) of the kinase, and could even favour it, binding to the high-affinity site is not. Stochastic simulations of a biochemical model showed that the existence of two such binding sites, one of them accessible only in the active, open conformation, would be sufficient to explain calmodulin trapping by CaMKII. We can explain the effect of CaMKII autophosphorylation at Thr286 on calmodulin trapping: It stabilises the active state and therefore makes the high-affinity binding site accessible. Crucially, a model with only one binding site where calmodulin binding and CaMKII inactivation are strictly mutually exclusive cannot reproduce calmodulin trapping. One of the predictions of our study is that calmodulin binding in itself is not sufficient for CaMKII activation, although high-affinity binding of calmodulin is

The Molecular Pages of the mesotelencephalic dopamine consortium (DopaNet)

BACKGROUND: DopaNet is a Systems Biology initiative that aims to investigate precisely and quantitatively all the aspects of neurotransmission in a specific neuronal system, the mesotelencephalic dopamine system. The project should lead to large-scale models of molecular and cellular processes involved in neuronal signaling. A prerequisite is the proper storage of knowledge coming from the literature. METHODS: DopaNet Molecular Pages are highly structured descriptions of quantitative parameters related to a specific molecular complex involved in neuronal signal processing. A Molecular Page is built by maintainers who are experts in the field, and responsible for the quality of the page content. Each piece of data is identified by a specific ontology code, annotated (method of acquisition, species, etc.) and linked to the relevant bibliography. The Molecular Pages are stored as XML files, and processed through the DopaNet Web Service, which provides functionalities to edit the Molecular Pages, to cross-link the Pages and generate the public display, and to search them. CONCLUSIONS: DopaNet Molecular Pages are one of the core resources of the DopaNet project but should be of widespread utility in the field of Systems Neurobiology

Modulation of calmodulin lobes by different targets: an allosteric model with hemiconcerted conformational transitions

Calmodulin, the ubiquitous calcium-activated second messenger in eukaryotes, is an extremely versatile molecule involved in many biological processes: muscular contraction, synaptic plasticity, circadian rhythm, and cell cycle, among others. The protein is structurally organised into two globular lobes, joined by a flexible linker. Calcium modulates calmodulin activity by favoring a conformational transition of each lobe from a closed conformation to an open conformation. Most targets have a strong preference for one conformation over the other, and depending on the free calcium concentration in a cell, particular sets of targets will preferentially interact with calmodulin. In turn, targets can increase or decrease the calcium affinity of the calmodulin molecules to which they bind. Interestingly, experiments with the tryptic fragments showed that most targets have a much lower affinity for the N-lobe than for the C-lobe. Hence, the latter predominates in the formation of most calmodulin-target complexes. We showed that a relatively simple allosteric mechanism, based the classic MWC model, can capture the observed modulation of both the isolated C-lobe, and intact calmodulin, by individual targets. Moreover, our model can be naturally extended to study how the calcium affinity of a single pool of calmodulin is modulated by a mixture of competing targets in vivo

Systems Biology Markup Language (SBML) Level 2: Structures and Facilities for Model Definitions

We present the Systems Biology Markup Language (SBML) SBML Level 2 Version 2, a model representation formalism for systems biology. SBML is oriented towards describing systems of biochemical reactions of the sort common in research on a number of topics, including cell signaling pathways, metabolic pathways, biochemical reactions, gene regulation, and many others. SBML is defined in a neutral fashion with respect to programming languages and software encoding; however, it is primarily oriented towards allowing models to be encoded using XML, the eXtensible Markup Language (Bosak and Bray, 1999; Bray et al., 2000). This document contains many examples of SBML models written in XML, as well as the text of an XML Schema (Biron and Malhotra, 2000; Fallside, 2000; Thompson et al., 2000) that defines SBML Level 2 Version 2. A downloadable copy of the XML Schema and other related documents and software are also available from the SBML project web site, http://sbml.org/. The SBML project is not an attempt to define a universal language for representing quantitative models. The rapidly evolving views of biological function, coupled with the vigorous rates at which new computational techniques and individual tools are being developed today, are incompatible with a one-size-fits-all idea of a universal language. A more realistic alternative is to acknowledge the diversity of approaches and methods being explored by different software tool developers, and seek a common intermediate format—a lingua franca—enabling communication of the most essential aspects of the models. The definition of the model description language presented here does not specify how programs should communicate or read/write SBML. We assume that for a simulation program to communicate a model encoded in SBML, the program will have to translate its internal data structures to and from SBML, use a suitable transmission medium and protocol, etc., but these issues are outside of the scope of this document

Ontologies for use in Systems Biology: SBO, KiSAO and TEDDY

The use of computational modelling in the description and analysis of biological systems is at the heart of Systems Biology. Besides the information stored in a core model, there is increasingly a need to provide additional semantic information: to identify model components, to assist in biological interpretation of models, to define simulation conditions and to describe simulation results. This information deficit can be addressed through the use of ontologies. We describe here three ontologies created specifically to address the needs of the Systems Biology community in each sub-division, and illustrate their practical use with the 'Repressilator' model (Elowitz and Leibler, 2000)

Neuronvisio: A Graphical User Interface with 3D Capabilities for NEURON

The NEURON simulation environment is a commonly used tool to perform electrical simulation of neurons and neuronal networks. The NEURON User Interface, based on the now discontinued InterViews library, provides some limited facilities to explore models and to plot their simulation results. Other limitations include the inability to generate a three-dimensional visualization, no standard mean to save the results of simulations, or to store the model geometry within the results. Neuronvisio (http://neuronvisio.org) aims to address these deficiencies through a set of well designed python APIs and provides an improved UI, allowing users to explore and interact with the model. Neuronvisio also facilitates access to previously published models, allowing users to browse, download, and locally run NEURON models stored in ModelDB. Neuronvisio uses the matplotlib library to plot simulation results and uses the HDF standard format to store simulation results. Neuronvisio can be viewed as an extension of NEURON, facilitating typical user workflows such as model browsing, selection, download, compilation, and simulation. The 3D viewer simplifies the exploration of complex model structure, while matplotlib permits the plotting of high-quality graphs. The newly introduced ability of saving numerical results allows users to perform additional analysis on their previous simulations

Ligand-dependent opening of the multiple AMPA receptor conductance states: a concerted model

Modulation of the properties of AMPA receptors at the post-synaptic membrane is one of the main suggested mechanisms behind synaptic plasticity in the central nervous system of vertebrates. Electrophysiological recordings of single channels stimulated with agonists showed that both recombinant and native AMPA receptors visit multiple conductance states in an agonist concentration dependent manner. We propose an allosteric model of the multiple conductance states based on concerted conformational transitions of the four subunits, as an iris diaphragm. Our model predicts that the thermodynamic behaviour of the conductance states upon full and partial agonist stimulations can be described with increased affinity of receptors as they progress to higher conductance states. The model also predicts existence of AMPA receptors in non-liganded conductive substates. However, spontaneous openings probability decreases with increasing conductances. Finally, we predict that the large conductance states are stabilized within the rise phase of a whole-cell EPSC in glutamatergic hippocampal neurons. Our model provides a mechanistic link between ligand concentration and conductance states that can explain thermodynamic and kinetic features of AMPA receptor gating.Comment: 4 figures, models available on demand. They will be published by BioModels Database upon publication of the articl